Effects of voids on the reconstruction of the equation of state of Dark Energy 
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We quantify the effects of the voids known to exist in the Universe upon the reconstruction of 
the dark energy equation of state w. We show that the effect can start to be comparable with some 
of the other errors taken into account when analysing supernova data, depending strongly upon 
the low redshift cut-off used in the sample. For the supernova data alone, the error induced in the 
reconstruction of w is much larger than the percent level. When the Baryonic Acoustic Oscillations 
and the Cosmic Microwave Background data are included in the fit, the effect of the voids upon the 
determination of w is much lessened, but is not much smaller than some of the other errors taken 
into consideration when performing such fits. We also look at the effect of voids upon the estimation 
of the equation of state when we allow w to vary over time and show that even when supernova, 
Cosmic Microwave Background and Baryonic Acoustic Oscillations data are used to constrain the 
equation of state, the best fit points in parameter space can change at the 10% level due to the 
presence of voids, and error-bars increase significantly. 



1. INTRODUCTION 

The current body of astronomical data makes it very 
difficult to live without dark energy. Combination of the 
supernova data with the 7 year results from WMAP and 
the Sloan Digital Sky Survey (SDSS) data on Baryonic 
Acoustic Oscillations combine almost seamlessly with ob- 
servations of the Hubble constant and the light element 
abundances to paint a picture which has become known 
as the standard cosmological model [lJ-Q- The Universe 
is not homogeneous however and observations of galaxy 
clustering, N-body simulations and the simplest theoret- 
ical considerations all predict that the majority of the 
volume of today's Universe is considerably underdense 
in terms of matter, which forms roughly shaped denser 
walls surrounding the underdense voids The Uni- 

verse is made up therefore mostly of the voids between 
overdense regions, a situation which has lead to evocative 
comparisons of the Universe with both a froth of bubbles 
and with Swiss Cheese [HI E| ■ 

Such underdense regions have been studied in detail 
to understand their effects on cosmological observables 
an d to investigate whether they can help explain 
the various observations of the Universe without dark 
energy - this might be possible if we are located 

close to the centre of a very large void [UUll- Such a 
void would typically have to be larger than 100 Mpc in 
radius, in comparison to the voids which have definitely 
been observed in the Universe which vary in size between 
very roughly 5-20 Mpc. In order to create such large 
voids, it would be necessary to have a very non-trivial 
power spectrum of perturbations f2~H . [25| . In this work 
we attempt to reach a much less ambitious goal and we 
only consider the voids which we have evidence to believe 
actually exist. 

The nature of dark energy is unknown. In particular, 
it is not known if its energy density remains completely 
constant over time or whether it is growing or decreasing. 
If dark energy can be modelled as a perfect fluid then any 
such evolution is determined by the equation of state of 



the fluid w which relates the energy density and the pres- 
sure, P = wp. By fitting supernova data, we are able to 
constrain that equation of state and future investigations 
hope to pin it down with increasing accuracy. 

Voids can be considered not only to provide additional 
constraints on w [HI, [27j , but also as potential sources of 
error. In this study we estimate what the effect of voids 
will be upon the reconstruction of the equation of state of 
dark energy, w, from the supernova data. In order to do 
this, we look at the evolution of voids in a Universe that 
contains collisionless matter and a cosmological constant. 
We consider distant voids in which supernovae explode 
and look at the effect upon photons as they climb out of 
their host void walls and move towards us. The presence 
of the voids increases the scatter in the luminosity red- 
shift plane - we aim to estimate this error and see how 
it compares with other errors on the reconstruction of w. 
We do not expect the error to be large since the fluctu- 
ations in the gravitational potential will be very small, 
less than a percent. However, since the velocity disper- 
sion of supernovae relative to the Hubble flow can create 
an observable effect, we feel it is worthwhile investigating 
the magnitude of the error due to these smaller voids. 

We do this in a way which is not entirely self-consistent 
- this shortcoming lies in the fact that we evolve voids in 
a ACDM Universe and then apply the estimated error 
to other kinds of dark energy where w ^= — 1. The rea- 
son we do this is that although it is in principle possible 
to study voids evolving in a background of dark energy 
with w 7^ — 1, the evolution of such voids involves the 
solution of a complex set of partial differential equations, 
whereas voids evolving in a cosmological constant back- 
ground are the result of ordinary differential equations 
and are therefore simpler to analyse numerically. We 
hope to estimate the magnitude of the error one expects 
from these walls in the knowledge that when one includes 
Baryonic Acoustic Oscillations (BAO) data and the Cos- 
mic Microwave Background (CMB), the equation of state 
w is constrained to lie quite close to — 1 anyway. We note 
that much more complicated analyses might be done in 
future if it becomes necessary. 
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In the next section we will outline the equations that 
we evolve as well as the functional form of the void models 
we choose. We will then describe our numerical proce- 
dures and present our results before making conclusions. 



2. THE LTB EQUATIONS AND VOIDS 

Although we will assume that all the matter in the Uni- 
verse has evolved into voids and void walls, our interest 
is focused on the hole nearest to the source only. We as- 
sume that the effects of complete holes situated between 
the source and the observer cancel out along the light 
path and we neglect the void that we are located in since 
all photons must arrive into the same current potential 
well. We construct our Swiss-Cheese model of the uni- 
verse by combining two components: the 'cheese', where 
the metric is described by a spatially flat Friedmann- 
Robertson- Walker (FRW) solution and the 'holes', where 
the matter is mainly underdense by volume and we use 
the Lemaitre-Tolman-Bondi solution. Unlike the idea 
that we are at the centre of a very large void which is 
used to live without dark energy, we will therefore be in 
the cheese, looking at distant holes. 



A. The Model: Lemaitre-Tolman-Bondi Universe 



Both of these functions are completely defined therefore 
by the choice of an initial density profile p(r, £ltb), where 
£ltb = ^ltb (r) refers to the beginning of the LTB evolu- 
tion and is set to a constant for simplicity. We then have 
to choose E(r) in order to match the flat FRW model at 
the boundary of the hole. This choice also ensures that 
the average density inside the hole equals the one out- 
side, so that an observer situated in the cheese would not 
be aware, locally, of the presence of the hole. 

We intend to test our model by sending photons 
through the voids and measuring redshifts. To do so, 
we must integrate the equation for their radial trajectory 



dt(r) 
dr 



R'(r,t(r)) 
(l + 2E(r)) 1/2 ' 



(5) 



To repeat and clarify, in contrast to other models based 
on the Swiss-Cheese universe, we assume that the ob- 
server, although occupying no particular position in space 
with respect to the holes, is always situated in the cheese; 
we then consider only the effects caused by the furthest 
hole where the supernovae are located (see Subsection 
3. A for more details) and hence the geodesies originate 
since in this study we want to investigate the effect of the 
void at the origin of the photons rather than when they 
arrive here. 



Our model is based on the Lemaitre-Tolman-Bondi 
(LTB) metric, a spherically symmetric solution of Ein- 
stein's equations, which can be written as (in units where 
c=l) 



ds 2 



-dt 2 + S 2 (r, t)dr 2 + R 2 (r, t) (d9 2 + sin 2 0d<p 2 ) ( 1 ) 



where we use comoving coordinates (r, 0, <p) and proper 
time t. For pressureless dust and a cosmological constant, 
Einstein's equations imply the following constraints: 



S 2 (r,t) 



1 . 2 GM{r) 
2 R 



R{r,t) 



R' 2 {r,t) 
1 + 2E(r) 



--AR 2 
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E(r), 



47rp(r, t) 



M'(r) 
R'(r,t)R 2 (r,t) : 



(2) 



(3) 



(4) 



where a dot stands for partial derivative with respect 
to t and a prime with respect to r; p(r, t) is the energy 
density of the matter, G is Newton's constant and A is the 
cosmological constant. To specify the model we intend 
to use, we have to define the two arbitary functions E(r), 
corresponding to the spatial curvature, and M(r), which 
is simply the mass integrated within a comoving radial 
coordinate r: 



M(r) 



4tt f p(r,t)R 2 R'dr. 
Jo 



B. Parameters and Initial Conditions 

To characterise the LTB model we intend to use, we 
have to define the initial density profile p(r, £ltb)- We 
then will build the two arbitrary functions E(r) and M(r) 
based on this profile. The expression we choose for it is 
based on Rostov's parametrisation [28[ where the density 
profile (Fig. [T]) is defined as follows: 



p(r,t ) = p(t Q )x 
{Ai + Aitanh [a (r 



r 1 )]-A 3 tanh\fi(r-r 2 )]}, (6) 



for r < r/j, where is the radius of the void; and 
p(r,to) = p(to) for r > Th- Given {r\,r?), the values 
of the coefficients {A\,Ai,A-$) and (a,/?) are chosen so 
that p{rh, to) = p(to) and the integrated mass inside the 
hole M(r h ) = ^ f h p(r,t )r 2 dr = f7rp(tb)rjj. 

The function E(r) plays the role of the spatial curva- 
ture in the LTB solution. It is defined using Eq. [3] and 
considering the initial time of the LTB evolution £ltb • 
At this early time, R(r, £ltb) = a(iLTB)^ = oltb?' and 
we have, in real units: 



E(r) 



LTB "LTB 



M(r) 



47T al TB rp(t LTB ) 



(7) 



where p(t) is the average energy density at time t. 

In order to guarantee that the chosen parametrisation 
of the voids is consistent with observational constraints, 
we compare the density profiles with observed values of 
over- and underdensities at different times of the LTB 
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Figure 1: Initial density profile define at t — £ltb and showing 
the shape of the cosmological voids and of its surrounding 
walls of denser matter. 
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Table I: Chosen parametrisation of the voids, taking into ac- 
count the constraints at initial time £ltb and final time t n ow 



evolution: at initial time t = £ltb and at final time 
t = t naw . We set tLTB to be the time of last scatter- 
ing, turB — t(z = 1100). So the constraint at initial time 
comes from the scale of temperature fluctuation in the 
Cosmic Microwave Background (CMB) AT/T « 10" 5 , 
which corresponds to a variation in density of the order 
of 10~ 4 . 

At the other extremity, when t = i now , the guideline 
for density profiles is given by observations of the present 
underdensities in the matter distribution of our universe. 
Recent studies (e.g. show that the voids that we 

see nowadays probably correspond to regions where the 
density is ~ 20% of the mean cosmic density. 

Further constraints on the size of the voids can be de- 
duced from the study of non-linear late-time integrated 
Sachs- Wolfe effects on the CMB [H|. 

Ensuring that these constraints are taken into account 
defines the spatial geometry of the voids. The chosen 
parameters are given in Table |TJ 



3. SIMULATION PROCESS AND RESULTS 
A. Obtaining Redshifts and Distance Moduli 

In order to compare the various Universes with dark 
energy possessing different equations of state, we need to 
build Hubble diagrams which allow us to analyse the dis- 
crepancies between models and Type la Supernova data. 
We therefore obtain the redshift of distant sources situ- 
ated around the void and then calculate the correspond- 
ing distance modulus. 

We evolve a single void up to some random time and 
then choose a random position in the void according to 
the weighting of matter density. We assume that the 
number density of supernovae explosions is proportional 
to the density and therefore more likely to occur in the 
void wall. We follow the path of a photon emerging from 
that void, with initial coordinates {r,i n ,ti n ), until the 
value of t(r) reaches the present time, t(ra na i) = t now . 
We then send a second photon, with initial coordinates 
( r 'in = r in,t' in = t ln + At in ), until the comoving distance 
it covered corresponds to rg na i. At this point, we can 
calculate At nna i = t'(Yfj na i) — < now and thus the redshift 
z(rj n ,ti n ) = Atfi na i/A<i n — 1. The luminosity distance is 
then given by 

d L = (1 + z) 2 i?(r fina i, t now ) = a r obs (l + z), 

where ao is the value of the scale factor today and 
''obs = '"fin — Hri is the distance between the source and 
the observer. We repeat this process, generating a cloud 
on the distance modulus / redshift plane which takes the 
form of a blurry line with scatter around the normal /i — z 
relationship in ACDM. 



B. Statistics: Binning Process 

The aim of the binning process is to gather series of 
points given by simulations in order to deduct sensible 
estimates of the error on the redshift and luminosity dis- 
tance from them. Given the redshift range covered by 
the supernova data, 0.01 < z < 2, we choose to dis- 
tribute the 10 5 data points between 200 bins of width 
Az = 0.01 (Fig. [2]). For each bin obtained, the mean 
and standard deviation are calculated. We thus end up 
with a new series of points and corresponding errorbars 
which we can now add to the existing known errorbars on 
the supernovae and see how the cosmological conclusions 
change. 

The Type la Supernovae data used in this study are 
extracted from the Union2 Compilation [l| of the Su- 
pernova Cosmology Project, which contains 557 sources 
drawn from 17 datasets. We base our comparison on the 
file called "Union2 Compilation Magnitude vs. Redshift 
Table" [33[, which contains supernova name, redshift, dis- 
tance modulus, and distance modulus error. To the latter 
we add the systematic error due to the peculiar veloci- 
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Figure 2: Illustration of the binning process allowing us to 
deal with the statistical nature of our results. The width of 
each bin is fixed in redshift, Az — 0.01. 

ties of the host galaxies, a v = 400 km/s, converted into 
a systematic error on the redshift of Sz = 0.00132. 

Once all data are binned and all errors known, we can 
compute x 2 values for every combination (fl m ,ft\). In 
parallel, we do the same for all FRW universes defined 
by the same density parameters. 



C. Comparison of Errors 

To compute consistent x 2 values, we ought to take into 
account all the errors and uncertainties adding up for 
each source considered and which can be split into two 
categories: error on distance modulus, a^, and binning 
process uncertainty, ctbp, so that 

2 x ^ [^nican(^) /^meas(^)] 
SNe ^ + ^BP 

Error on Distance Modulus | This first error gathers 
the variety of measurement errors due to the nature of the 
sources and the way their luminosity is evaluated. Ac- 
cordingly to the Union2 original paper [l[ , we sum in the 
component errors coming from the following causes: light 
curve fitting (through covariant matrix), galaxy peculiar 
velocities, Galactic extinction, gravitational lensing and a 
floating dispersion term (for potential sample-dependent 
systematic errors). 

All these errors are included in the error on distance mod- 
ulus given by the Union2 Compilation data. 

Binning Process Uncertainty | The second group of 
errors relate to the spread induced by the voids them- 
selves. 

Here the error comes from the binning process itself, 
where we combine results of simulation for a large num- 
ber of different sets of initial conditions. For each bin 



0,05 0.1 0.15 0.2 025 0.3 035 0.4 0.45 0.5 

Figure 3: Qm vs dark energy equation of state w when fitting 
the Union2 data set only assuming a flat Universe. The shad- 
ing represents probability from the % 2 cumulative distribution 
function. 

created at a given redshift z (see Subsection 3.B), we 
determine the standard deviation of the data by evaluat- 
ing the spans from the mean value that include 68% of 
all the points in the bin. The corresponding deviation, 
cbp(z), is then simply added in quadrature to the total 
uncertainty. 



D. Effect on the estimates of the Equation of State w 

Having estimated the errors due to the growth of voids 
in the universe, we are in a position to find out how this 
effects the reconstruction of the equation of state. To do 
this we perform a Monte Carlo Markov Chain (MCMC) 
fitting to the data. We assume a flat Universe and gener- 
ate values of (and therefore Qde = 1 — ^m) an d W 
randomly. For all fits in this paper we remove uncertainty 
in the Hubble constant by choosing the best fit value for 
each value of Hm- We also calculate the probability from 
the x 2 values using the % 2 cumulative distribution func- 
tion rather than relative x 2 values as an estimation of 
maximum likelihood. 

First we fit the supernova data alone. Typical favoured 
regions in the J7m — w plane can be seen in figure [3] 
Table [XT] shows the projected constraints on the equa- 
tion of state with and without the presence of the ad- 
ditional error due to voids. We have done MCMC runs 
for the supernova data sets with various minimum red- 
shifts. When fitting to data, it is normal to cut out very 
small redshift supernova since they are the data points 
most vulnerable to large redshift errors such as peculiar 
velocities. Given that we expect the void effect to be 
at most the same magnitude as these other errors, we 
perform fits to the data with and without the voids ne- 
glecting supernovae below different minimum redshifts 
- z mm = 0.01,0.02,0.03 and 0.04. As expected, the dis- 
crepancy between the simulations with and without voids 
decreases as we increase the minimum redshift. 
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Table II: The best fit for the equation of state of dark energy w 
to the Union2 supernova data set alone with the 67% errors. 
The fits are done for different minimum redshifts and with 
and without the additional errors created by the presence of 
voids. The fourth Column is the difference between the two 
best fit values. Note the general trend of the difference in 
errors going down as we cut out lower redshift supernovae 
(apart from when we go from z min = 0.02 to z min = 0.03 
which must be a statistical fluctuation). 
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Table III: The best fit for the equation of state of dark en- 
ergy w to the Union2 supernova data set, CMB and BAO, 
as described in the text, with the 67% errors. The fits are 
done for different minimum redshifts and with and without 
the additional errors created by the presence of voids. 
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Table IV: The best fit for the dynamic equation of state of 
dark energy w = wo + w a yrj to the Union2 supernova data 
set, CMB and BAO, with the 67% errors. The fits are done 
for different minimum redshifts and with and without the 
additional errors created by the presence of voids. 



Figure 4: Qm vs dark energy equation of state w when fitting 
the Union2 data set, Baryonic Acoustic Oscillations and the 
shift parameter and acoustic scale of the CMB, again assum- 
ing a flat Universe. The shading represents probability from 
the x 2 cumulative distribution function. 



The constraint on the value of w when only including 
the supernova data is very weak and best fit values are 
quite far from w = — 1. This regime is therefore far 
from the region of validity of our analysis - the growth of 
voids in such a Universe might be quite different to the 
LTB case due to structure formation in the dark energy 
and anisotropic pressure which in general is a challenge 
to calculate. However, we make no apologies for this as 
we only seek to highlight the possible magnitude of void 
effects upon the reconstruction of the equation of state 
in this work. 

Next we include the effective distance measure Dy 
which is obtained from the SDSS Baryonic Acoustic Os- 
cillations data [!, [3(| and the WMAP constraint on dark 
energy from the shift parameter R and the acoustic scale 
I A |2|. The favoured regions on the f2 — w plane can be 
seen in figure @] These lead to much tighter constraints 
on possible values of w and the difference between the 
cases with and without voids is much smaller (Table IITT]) . 

We note here that we obtain slightly different con- 



straints on the equation of state to the analyses carried 
out in [l|, which is probably mostly due to the way we 
marginalise over the Hubble constant. The basic pre- 
ferred regions are compatible in the two studies. 

E. Effect on the estimates of the dynamic Equation 
of State w(z) 

So far, we have considered a constant value for the 
equation of state parameter w. ft is however possi- 
ble that the dark energy equation of state might vary 
over time and therefore redshift. A very commonly used 
parametrization of dynamic dark energy is 

z 

w = w + w a — , 

l + z 

based on works by Chevallier, Polarski and Linder [3ll . 
[32j]. Following the same procedure as before, we fit the 
supernovae data while taking into account the constraints 
given by the CMB and the BAO. The results for the best 
fits are given in Table H*Vl 

Even with the extra degree of freedom, the fits ob- 
tained when including the whole set of supernovae down 
to arbitrary small redshifts are not fantastic, the prob- 
ability obtained from the x 2 distribution of the best fit 
point in parameter space being less than 50%. However, 



(i 



the situation changes dramatically depending on the red- 
shift considered for the cut-off: the best fit probability 
comes close to 1 with z m i n = 0.04. Fig. [5] shows the 
contour plots obtained for the different values of z m i n 
studied here. 

It is clear again that if we include the lowest redshift 
supernovae, the error which is introduced due to the pres- 
ence of voids has a considerbale impact on the best-fitting 
regions in w a — w a parameter space. 



4. COMMENTS AND CONCLUSIONS 
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There are more than 100 supernovae below z = 0.04 in 
the Union2 data set so it is not surprising that including 
or not including them can have a considerable effect on 
the reconstruction of the equation of state of dark energy 
w. In the same way, changing their error bars changes 
their pull on the best fit values of w. 

In this work, we have shown that the growth of voids in 
the Universe can have a small effect upon the reconstruc- 
tion of the equation of state of dark energy by inducing 
scatter into the redshift-luminosity data plane. If one 
considers only the supernova data, and we try to con- 
strain models of dark energy with a constant equation of 
state, the effect of the voids may be relatively large, up 
to around 10% of the central value of w. becoming less as 
one increases the minimum redshift supernovae included 
in the fit. 

When one includes the constraints from the CMB and 
from Baryonic Acoustic Oscillations, the magnitude of 
the effect is much smaller, below the percentage level. 
However, the induced error does start to compare with 
some of the other smaller errors that are considered by 
the Union2 team which are listed in table 9 of [l| . This er- 
ror is not important now and may turn out to be dwarfed 
by other systematics even in future studies using data 
from new observational programs. Conversely, one can 
envisage a situation where it would have an important ef- 
fect upon the understanding of dark energy and we would 
need to parametrise and quantify our uncertainty. 

If we consider dark energy with an equation of state 
which is allowed to vary, i.e. w(z) = w a + w a z/(l + 
z), we find that the effect of the uncertainty introduced 
due to the voids is amplified and the best fit values of 
w and w a can change around 10%, unless we neglect 
all supernovae at redshifts lower than z ~ 0.04. The 
magnitude of the errors around these best fit regions also 
increase significantly. 

Another approach to understand the distribution of 
these voids and their effect is to try and map out the 
local voids using peculiar velocity flows in an attempt to 
quantify their effect. This seems difficult : a redshift of 
z ~ 0.04 corresponds to a comoving distance of roughly 
150 Mpc. The size of typical voids in the Universe today 
is considerably smaller than this, of the order of 10-20 
Mpc in radius. We should therefore expect many hun- 
dreds of voids to be present in this low redshift part of the 
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Figure 5: Dark energy equation of state parameters wo versus 
w a , with z m in = 0.01,0.02,0.03 and 0.04 respectively, when 
fitting the Union2 data set, Baryonic Acoustic Oscillations 
and the shift parameter and acoustic scale of the CMB, always 
assuming a flat Universe. The contours correspond to the 
probabilities listed in the individual graph titles. 
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Hubble diagram and it seems impossible to trace them 
given our lack of understanding of bias and the history 
of structure formation in the local Universe. What we do 
know is that there are voids, their gravitational potential 
may have implications for the reconstruction of w and we 
can in principle quantify this uncertainty. 
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